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Technical  Objectives 

The  objectives  of  the  project  are  three-fold:  (1)  Development  of  first-principles  simulator  to  understand 
the  fundamental  physical  processes  underlying  the  performance  of  the  solar  cells.  (2)  Identifying  key 
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Technical  Approach 
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in  the  simulations.  More  specifically,  this  approach  is  based  on  the  time-domain  ab  initio  non-adiabatic 
molecular  dynamics  for  simulating  phonon-assisted  electron  transitions.  The  electronic  energy  levels  and 
their  transition  rates  are  determined  from  the  density  functional  theory  (DFT),  taking  into  account  of  both 
inter-molecule  and  intra-molecule  contributions  and  treating  dynamic  and  static  disorder  at  an  equal 
footing.  This  multiscale  framework  incorporates  quantum  mechanical  determination  of  electron  transition 
rates,  mesoscale  master  equations  for  carrier  mobility  and  continuum  models  for  current-voltage 
characteristics  tailored  for  experimental  measurements.  For  interfacial  exciton  dissociation,  we  have 
developed  a  first-principle  approach  based  on  the  time-dependent  density  functional  theory  (TDDFT)  to 
describe  exciton  states,  including  energy  levels  and  many-body  wave  functions.  The  non-adiabatic  ab 
initio  molecular  dynamics  (MD)  is  used  to  determine  the  phonon-assisted  transition  rates  between 
localized  exciton  states.  In  conjunction  with  the  Monte  Carlo  and  surface-hopping  methods,  this  approach 
can  simulate  exciton  dynamics  in  organic  semiconductors  at  different  temperatures. 
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Figure  1:  (a)  Charge  separation  efficiency  as  a  function  of  L  with  M  =  40.  (b)  Charge  separation  efficiency  as  a  function  of  M 
with  L  =  41.  Hole  (c)  and  electron  (d)  diffusion  distance  along  x,  y  and  z  axes  averaged  over  all  KMC  trajectories  with  energetic 
disorder  of  10  meV.  The  diffusion  distance  for  the  hole  (e)  and  electron  (f)  along  x,  y  and  z  axes  averaged  over  all  KMC 
trajectories  with  energetic  disorder  of  70  meV. 

Energy  independent,  yet  higher  than  90%  internal  quantum  efficiency  (IQE)  has  been  observed  in  many 
organic  photovoltaics  (OPVs).  However,  its  physical  origin  remains  largely  unknown  and  controversial. 
Hypothesis  that  the  lowest  charge-transfer  (CT)  state  may  be  weakly  bound  at  the  interface  has  been 
proposed  to  rationalize  the  experimental  observations.  In  this  work,  we  study  the  nature  of  the  lowest- 
energy  CT  (CTi)  state,  and  show  conclusively  that  CTi  state  is  localized  in  typical  OPVs.  The  electronic 
couplings  in  the  donor  and  acceptor  are  found  to  determine  the  localization  of  CTi  state.  We  examine 
geminate  recombination  of  CTi  state  and  estimate  its  lifetime  from  first-principles.  We  identify  the 
vibrational  modes  that  contribute  to  the  geminate  recombination.  With  materials  parameters  determined 
from  first-principles  and  experiments,  we  carry  out  Kinetic  Monte  Carlo  simulations  to  examine  charge 
separation  of  localized  CTi  state.  We  find  that  the  localized  CTi  state  can  indeed  yield  efficient  charge 
separation  with  IQE  higher  than  90%.  Dynamic  disorder  and  configuration  entropy  can  provide  energetic 
and  entropy  driving  force  for  the  charge  separation.  Charge  separation  efficiency  depends  more 
sensitively  on  the  dimension  and  crystallinity  of  the  acceptor  parallel  to  the  interface  than  normal  to  the 
interface.  Reorganization  energy  is  found  to  be  the  most  important  materials  parameter  for  charge 
separation  and  lowering  the  reorganization  energy  of  the  donor  should  be  pursued  in  materials  design. 

2.  Interfacial  Oxygen  Vacancies  as  a  Potential  Cause  of  Hysteresis  in  Perovskite  Solar  Cells 
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Figure  2:  (a)  Electron  transfer  fraction  (%)  from  the  perovskite  to  the  TiCE  conduction  bands  as  a  function  of  time  for  different 
densities  of  oxygen  vacancy  (VO)  in  TiCE.  (b)  Electron  transfer  as  a  function  of  time  during  the  recombination  process  for 
different  densities  of  VO  in  TiC>2.  (c)  Electron  transfer  fraction  from  the  perovskite  to  the  TiCE  conduction  bands  as  a  function  of 
time  for  intact  Ti02  and  Ti02  with  VO  at  different  sites,  (d)  Electron  transfer  as  a  function  of  time  during  the  recombination 
process  for  intact  Ti02  and  Ti02  with  different  VO  sites. 

Organometal  halide  perovskite  solar  cells  (PSCs)  have  emerged  as  one  of  the  most  promising 
photovoltaic  technologies  with  efficiencies  exceeding  20.3%.  However,  device  stability  problems 
including  hysteresis  in  current-voltage  scans  must  be  resolved  before  the  commercialization  of  PSCs. 
Transient  absorption  measurements  and  first-principles  calculations  indicate  that  the  migration  of  oxygen 
vacancies  in  the  TiC>2  electrode  under  electric  field  during  voltage  scans  contributes  to  the  anomalous 
hysteresis  in  PSCs.  The  accumulation  of  oxygen  vacancies  at  the  electrode/perovskite  interface  slows 
down  charge  extraction  while  significantly  speeding  up  charge  recombination  at  the  interface.  Moreover, 
first-principles  non-adiabatic  molecular  dynamics  simulations  reveal  that  the  charge  recombination  rates 
at  the  interface  depend  sensitively  (with  1  order  of  magnitude  difference)  on  the  locations  of  oxygen 
vacancies.  By  intentionally  reducing  oxygen  vacancies  in  the  TiCB  electrode,  we  substantially  suppress 
unfavorable  hysteresis  in  the  PSC  devices.  This  work  establishes  a  firm  link  between  microscopic 
interfacial  structure  and  macroscopic  device  performance  of  PSCs,  providing  important  clues  for  future 
device  design  and  optimization. 


3.  Charge  Stripe  Formation  in  Molecular 
Charge  Separation 


Ferroelectric  Organohalide  Perovskites  for  Efficient 


Figure  3:  Charge  density  of  the  lowest  energy  exciton  in  MAPbfi  (left)  and  NaPbfi  (right)  viewed  along  the  z  direction.  The 
electron  and  hole  density  is  colored  in  blue  and  yellow,  respectively. 
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Despite  rapid  progress  in  the  efficiency  of  organohalide  perovskite  based  solar  cells,  physical  mechanisms 
underlying  their  efficient  charge  separation  and  slow  charge  recombination  still  elude  us.  Here  we  provide 
direct  evidence  of  spontaneous  charge  separation  via  first-principles  simulations.  The  excitons  are 
predicted  to  self-organize  into  stripes  of  photo-excited  electrons  and  holes,  spatially  separated  as  effective 
channels  for  charge  transport.  The  rotation  of  organic  cations  deforms  the  inorganic  framework,  and  as 
the  deformation  reaches  a  critical  value,  a  direct  band  gap  transforms  to  an  indirect  one,  and  the  photo- 
excited  electrons  rotate  in  alignment  with  the  deformation-induced  electric  fields.  The  latter  triggers  Stark 
effect  which  in  turn  leads  to  the  formation  of  charge  stripes.  The  interplay  between  dynamic  disorder, 
ionic  bonding  and  polarization  is  responsible  for  the  formation  of  the  charge  stripes  and  the  indirect  band 
gap,  both  of  which  could  lead  to  efficient  charge  separation  and  reduced  charge  recombination  in  the 
organohalide  perovskites. 
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